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In this work we solve two new matrix models, using standard and new techniques. 

The two models are based on matrix ensembles not previously considered, namely 

the Dili generator ensemble. It is shown that, in the double scaling limit, their free 

C^ I energy has the same behavior as previous models describing oriented and unoriented 

surfaces. We also found an additional solution for one of the models. 
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Introduction 

Recent progress has been made in non-perturbative string theory via matrix models, where 
one represents triangulated random surfaces by large-A^ matrix integrals. By taking the 
so-called double scaling limit, one gets some insight in the physics of the model, associated 
with models of 2-D gravity. In this paper we look at two new such matrix models. 

The reasons for studying them are that they were part of a classification scheme, but 
also they provided an opportunity to extend the techniques for solving matrix models, 
namely via skew-orthogonal polynomials |]l| . We also had to solve them to compare with 
previous models to see if the physics revealed by the free energy is the same. Indeed, from 
its series expansion, we can say if the model describes unoriented surfaces in addition to 
oriented ones, and, by calculating an appropriate ratio, we can say if it is the same as 
previous models. 

Section one is devoted to a classification scheme for matrix models within the context 
of symmetric spaces. In section two and three, we show the calculations of two new models 
in D = 0. Finally, we conclude with a discussion and a comparison of the two models' 
free energies, and a comparison with other models. 

1 Classification of matrix models 

We will now introduce a classification of matrix models based on symmetric spaces. This is 
an interesting class of single-matrix models for which the reduction of the matrix integral 
allows for a solution by a polynomial method. 

Symmetric spaces are cosets of the form G/K. Technically speaking, those that we 



are considering are the simply connected Riemannian globally symmetric spaces P, p|, p 
The classical types are listed in table |1]. The two matrix models, in sections two and three, 
are based on the Dili generator ensemble. They correspond, for section two, to n odd, 
and for section three, to n even. 

There are a few ways to reduce a matrix integral, depending on the type of matrix 
or the matrix ensemble. For example, for the Hermitian matrix model, one can simply 
diagonalize the matrix as it is done in Metha's book [^]. The result is that we get an 
integral over the eigenvalues with a Jacobian in terms of these eigenvalues. The idea 
behind this reduction is to simplify the matrix integral. 

The symmetric spaces introduced in the previous subsection give us a systematic 



approach to identify the matrix ensemble for which a similar reduction is possible. For 
these spaces, however, the procedure of reducing the matrix integral is more complex. An 
analogy can be used to understand how the Jacobian and the set of integration parameters 
arise. Consider the volume element in flat space written in cartesian coordinates and in 
spherical coordinates, Ylf^idxi = dujdrr^^^. Spherical symmetry would allow us to 
integrate out the angular part du with r^~^ being the Jacobian and r the integration 
variable. For a symmetric space (matrix ensemble) the Haar measure can be decomposed 
in the same way, with the symmetry provided by the associated group. In doing so we get 
the Jacobian and also a set of parameters (which are not necessarily the eigenvalues of the 
matrices). The example above can be regarded as the measure for the generator ensemble 
of the coset SO{N + 1)/S0{N). The parameter of interest is r and du is the measure on 
SO{N), which yields a "trivial" integral if the original integrand is spherically symmetric. 
Such vector models were studied in [^]. For the technical details the reader is referred to 
Helgason's book 0. Let us just say here that there is a relation between the multiplicities 
of the roots and the exponents in the Jacobian. Table ^ gives some symmetric spaces 
and their associated Jacobian. Many of these models have been previously considered, 
including A I |§|, 0, corresponding to real symmetric matrices, A II [0, |[, corresponding 
to quaternionic real self dual matrices, as well as A III , BD I, and C II P, . 

It should be noted here that matrix models with Hermitian, antisymmetric, and 
quaternionic real self dual matrices, are not part of this classification. In fact, they 
are generator ensembles for the classical groups U{n),SO{n), Sp{2n), respectively. 

In the present paper we consider the D III generator ensemble, so only C I remains to 
be studied. 

2 First model — n odd 

We will first study the double-scaling limit of the D III generator ensemble with n odd 
(see table ^. The matrix integral always takes the form, 

Z = J dMe-f^'''^^\ (1) 

The Jacobian for this model is 

TV N 

:^=ni^'in(^'-^?)', (2) 

i i<j 



and we are going to make the analysis with the potential V = ax^ + h/2 x'^. The partition 
function can then be written as follows, 
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where we made the change of variables x"^ = y. For convenience of evaluation we will 
rewrite it as a determinant [H|. We have, 
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-PV{y) 



The P's are the usual orthogonal polynomials on the half-line, defined as 

Pk{y) = y'' + l.o. Po{y) = l 

POO 

/ dfiPk{y)Pi{y) = hk5i^k d^i = dye- 
Jo 

and from which we derive the following recursion relation, 

yPk = Pk+l + SkPk + RkPk-l, 



(5) 
(6) 

(7) 



with Rk = hk/hk-i. 

In order to solve for the free energy, we will need to know, first, the solutions for R 
and S. By considering recursion relations for dyPk and ydyPk, one finds a set of two 
self-consistent equations in terms of R, S, k, a, b, and (3 pUj] , 

4:k + l 



/3 



0^S2k + b{R2k + R2k+1 + "S*. 



2k) 



and 



(a + bS2k)ia + bS2k-i)R2k = j[S2k{a + bS2k) + b{R2k+i - R2k) 



P' 



(9) 



We used 2k indices because in the final recursion relation, we will need to know the 

solutions of R2N+1 and S2n+i- In the planar approximation, and choosing Ac = {k/ [3)c = 

1/4 and Re = I (together with the criticality condition dsX = ), we find all the critical 

values, 

Ac = 1/4 R, = l S, = 2 
a = l 6 =-1/6. 



(10) 



We choose the scahng solutions to be, 

S2N+1 = 2(i-r'^exp(-^r'^^)/) (11) 

R2N+1 = l-f3-'expi-^-r''^^)i9o + p-'9i + 92p-'' + ...) (12) 

where t is defined by t = (/5/4 — N)P~'^. We write A^ instead of k because we are only 
interested in the large N behavior of R and S (and so we replace A; by A^ in our recursion 
relations as well). In the R ansatz we added some more degrees of freedom by expanding 
the g function in powers of [3'". This will be necessary to get consistent solutions beyond 
leading order in jS'". Doing the same thing for / would only yield redundant equations. 

After inserting in the recursion relations, and working out the lowest order equations, 
we see that a consistent solution requires /i = 2z/. For the two recursion relations, the 
results are given, order by order, in table ^. 

From eq.(§), we obtain the values of the exponents, u = 1/5, and fi = 2/5, which is 
consistent with our previous relation between /i and z/. Also, using the relations between 
the g coefficients and /, we get, at order j3~'^~^'^ in (^, a differential equation for /, 

4t = /2-l/48/" (13) 

which is the well-known Painleve I equation, upon renormalization of t. At order /3~^ and 
j3~^~'^ , we simply get equalities (e.g. = 0), which come from our criticality requirement 
( first derivative of A with respect to S is zero). 

To solve the model, one must still solve for the double-scaling limit of the entire 
partition function. For this, it is clear that we have to find a recursion relation for 
y'^dyPniy) (as appears in the determinant H). After some algebraic manipulations, one 
finds 

y^dyPn = nPn+l + I ^Cn,n - Snj Pn + {PCn-l,n -U- l)i?„Pn„l + 

(3Cn-2,nRnRn-lPn~2 + PCn~3,nRnRn-lRn-2Pn-3, (14) 

where, for simplicity we wrote C^^i = J dfiPkPiy'^V'{y). With this recursion relation, the 
partition function is found to be, 

Z^+l = iN+l){< P^j^,P2N+l >-< P2N+1^P2N >)Zn + N{N+1) 
< P2N-2^P2N+1 > (< P2N~1^P2N > — < P2N,P2N-1 >)^Af-l 
-N{N + 1){N-1)< P^j,_,, P2N+I >< P^N^s, P2N > 

< P^j^_^, P2N-1 > Zm-2 - (A^ + 1) < P^^„i, P2N+1 > Yn, (15) 



where P' = y^dyP. Yn is an auxilliary partition function with a determinant similar to 
Zjv but where the last and third-to-last columns were removed. It was introduced to 
avoid an infinite number of terms in the previous relation, and satisfies 

Yi^ = N < P2N~2^ P'iN > 2^N-1 — N < -P2Af-35 ^^N > ^N-l- (16) 

We now define, for later convenience (we want to have a smooth planar limit as A^ ^ oo ), 
the following two ansatz, 

^N — ^ ,.^,, Aat - ^ — — , {l() 

ZN^iN(3hh2N-i ZN-iNf3bh2N-i 

and rewrite eqs.(|l^) and (|16D in terms of polynomials in W and X. After some algebraic 
manipulations, and expanding the brackets, we get 

WnWn-i 

Wn+iWnWn-I 7 {ci{S2N+l + 'S'2Ar) + h{R2N + R2N+I + R2N+2 + 'S'2^4.1 

~^^2N + S2N+1S2N) 7^ 7^) H ; (o. + b{S2N+l + S2N + S2N~l)) 

p p 

— Wiy^lR2N{R2N+l + R2N + R2N-I + 'S'ajv-l + 'S'2^ + S2N-1S2N) 

0'iS2N + 'S'2Af-l) -^ + R2nR2N-iR2N~2 = (18) 



b \ p J 

Xn+iWn ; {a + b{S2N+2 + S2N+1 + S2n)) + R2N+2XN = 0. (19) 

b 

In the planar limit, and using the previously found critical values for i?, S, etc, we find 
that the critical values are Wc = —1 (triple root of ITB|), while Xc remains undetermined. 



X it will be determined below when solving the recursion relations. 
To solve them, we choose the following ansatz, 

W^+i = -l + p-Pexpi-ip-''^)iho + p-''h+p-'"'h2 + ...) (20) 

at 

X^+i = X,-p-''exp{-ip-'^^){ko + p-'h+p-'''k2 + ...). (21) 

After insertion in eqs. ([TSD and (0) (at P~'^'^ order ), and using = /i = 2/5, a = p = 
V = 1/5, we find the results given in table §. We see that a consistent solution requires 
that Xc = 0. In doing so, we get two differential equations defining h^ and ko, 

= 6f + hoko + k'o (22) 

= -6/' - 12fho + hl + 3hoh'Q + h'^ - 6fko. (23) 



We used the fact that go = 2/ (table |^). Using the known solution for /, 

/ = 2ti/2 - — 1~2 _ ^^t-9/2 _ _ _ (24) 

•' 384 589824 ' ^ ' 

and power series solution for Hq and kQ,we find a set of algebraic equations that we can 

solve. We finally end up with four solutions, 

ho = ±2V3f/' + It-' T^r'^' + ... (25) 

ko = T2y3tV^ + V ± l^t-^/^ + . . . 



or 

ho 

ko 



±2V3t'/^ - Ir' T ^t-9/^ + . . . (26) 

T2V3t'/' + y ± ^t-^/^ + . . . . 



It turns out that we will only need ho in the solution of the free energy. Using the following 
relation for large N, 

exp {-d^F) = ^2 = R2N+iR2N-r7^, (27) 

the second derivative of the free energy is found to be 

F" ^ 4/ - h',. (28) 

Using the solution for / and the solution (EBf) for ho, we find, 



2 96 ^ ' 

In order to compare this result with other matrix models, we need a quantity independent 
of the scale of t and F". Such a quantity is the product of the first and the third coefficient 
in the series for F" divided by the square of the second coefficient. Here we find —13/9, 
which is the same result as in [^ |^ although the matrix ensemble in these papers is 
different. From previous solutions this is the expected ratio for the free energy of pure 



2D quantum gravity with oriented and unoriented surfaces. Doing the same with (pSf) , 
we get, 

F" = 8ti/2 ^ ^^-3/4 83^_2 

2 96 ^ ^ 

In that case, the universal ratio CoC2/cf yields 83/9, which differs from any known matrix 

models. One surprising result is that our two solutions differ only by t"^ (verified up to 

15th order of F" ), which is the term corresponding to the torus and the Klein bottle. 



3 Second model — n even 



We start again from the matrix integral ([^). But now, the Jacobian takes the form, 

N N 

J = Il\^^\E{x^-x',r. (31) 

i i<j 

We are doing the analysis, as usual, with the potential V = ax^ + h/2x^. The partition 
function can then be written as follows, 



N N 



■^0 i i<j 



i -f3(ay,+b/2y^) 



(32) 



We used, again, the substitution x^ = y. The complicated part is the Jacobian. For 
convenience of evaluation we rewrite it as a determinant [|T|. We have. 



N 






Qoivi) Qiivi) 

dyQoivi) dyQiiyi) 

Qo{y2) Qi{y2) 



Q2N-i{yi) 

dyQ2N-l{yi) 
Q2N-~l{y2) 



dyQoiyN) dyQi{yN) ■■■ dyQ2N-i{yN) 
where the Q's are Metha's skew-orthogonal polynomials, 

Qi{y) =y" + l.o., 



(33) 



and 



with 



< Q,,Q, >Q^^Jdye-^''^y\Q,Q'^ - Q\Q,) = g[,/2]%. 



Z2i,2i+1 — 1 
Z2i+l,2i = —1, 

all others being 0. With these definitions we can easily evaluate Z^r, 

N-l 

Zj, = 2-^iV! n q^■ 

1=0 



(34) 
(35) 

(36) 

(37) 



(3^ 



Here we used skew-orthogonal polynomials because with the orthogonal ones, Z^ 
cannot be evaluated (i.e. gives an infinite recursion relation for dP ). On the other hand. 
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with Q polynomials, Z^ is easy to find, but the problem is to establish recursion relations 
for them and the g's. Indeed, the only know recursion relation is infinite [|l|. The approach 
that we will follow is to relate the Q polynomials with the P polynomials (these are the 
usual orthogonal polynomials for which recursion relations are well-known, or at least, 
easy to find). 

We start with the general expansion of the P's in term of the Q's, 

P2i = Q2i + ^ilQli-l + ^i2Q2i-2 + • • • (39) 

P2i+1 = Q2i+1 + iilQ2i + C,i2Q2i-l + • ■ ■ • (40) 

In order to find odd and even ^'s and u 's, we consider four Q products. In the calculation 
we take the Q product of a P polynomial and a Q polynomial. After using the expansion 
of P's in terms of Q's, and rewriting the Q product in a P product (orthogonal polynomial 
relation), we get the following relations, 

Ci,2jqi-j = -< f3V'Q2i-2j,P2i+l>P --Q2i-2MP2^+liO) (41) 

U;i,2j-m-j = -<f3V'Q2i-2j,P2i>p--Q2i^2A0)P2i{0) (42) 

—^i,2jQi-j = -^ < PV Q2i-2j+l,P2i >P ~-^Q2i-2j+l{0)P2i{0) 

-^(2t + l)5joh2i (43) 

-^^,2j+lq^-j = -< /3V'Q2i-2j+l,P2i+l>P --Q2i-2j+l{0)P2i+l{0). (44) 

Considering the above equations for specific values of j , and defining the following quan- 
tities which have a smooth planar limit. 



W. -J^ Y - ^2^(Q) V. - ^2^+1 (0) 

' f3h2. ' P2.(o) '' P2m 

Pfc+i(0) , PkjO)' 
'~ P.(0) '"'- (3h ' 

we finally get, 

Wj ^ b 1 XiA2i+i 

R2i+1 2 2 Z2i 

^i2Wi_, 1 Xi_^A2i+l 



(45) 



(46) 



R2i+lR2iR2i-l 2 Z2iZ2i-lZ2i-2 

= a + hS2i - ojiih - XiA2i (4. 

uJiiWj-i ^ 1 Xi^iA2i 

R2iR2i-l 2 Z2i-lZ2i-2 



(47) 
(48) 
(49) 



-2Wi = hR2i+i - ^ii{a + hS2i) + [iiiUJii - i,2)h - Y^A^i - ^^'^ ^^ (50) 
-^isW^^-i ^ & 1 Yi_iA2i 

-R2J-R2J-I 2 2Z2j_l2'2j_2 



a + hS2i+i - Kii . (52) 

-K2J+I ^2i 

In the planar limit, and using critical values found in the previous section, we find, 
b -^ -1/6 \ = i/l3^llA S ^2 Z ^ -1 cj,,„ -^ uj^ 

a^l R^l A ^2/3 iln^in .„^ 

X ^ 1/2 W -^ 1/12 uJi^-2 6^2 ^1 ^ -3/2 ^ ^ 

F -> 1/2 UJ2 -^ 3, 

where the values for A and Z can be easily found using their explicit form. 

We can take a look at the partition function to see exactly what quantities we have 



to know. From (^), and using (^5]) , we get 



However, the ratio of Z's is related to F", so all we have to know is the ratio Wn IWn-i, 
something that we can easily find with a recursion relation for the 1^'s (of course we also 
have to find the differential equation satisfied by the function used in the W ansatz). 

To find this recursion relation we only need eqs.(PB|), (^), and (^5|). After some 
algebraic manipulations we get a recursion relation for Wn-, 

(- ^^^ "^ ^^ + hR2N+l){-^ + hR2N)[2WN - hR2N+l]WN-l 

2N y^ 

+R2N+i{a + bS2N){a + bS2N+i)[{-^WN~i) + -R2nR2N^i] = 0. (55) 

Using some information from table ^, which applied for recursion relations of P's, and 
using (j) = ^ = 2/5, p = u = 1/5, we find the coefficients listed in table |^. Replacing 
solutions for R2N+1, and W2N+1 in (^41), we get, 

F" = 4/ - h',. (56) 

Using a power series expansion for / and replacing in 6/ = /iq + /iq, we find 

ho = ±2V3t'/^ - -r' T ^r^/^ + . . . . (57) 

" 8 384 ^ ' 

Finally, we end up with the following solution for the second derivative of the free energy, 
F" = 4/ - h', ^ 8tV2 ^ %^,-3/4 _ 13^-2 ^ ^ ^ _ (58) 
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The appearance of a term 0{t~^'^) reveals that this model describes unoriented surfaces 
for sure, as well as oriented ones. To compare with previous model, we consider the ratio 
of coefficients. 



4 (^^^)2 9 ■ 

This is the same ratio as what was found in the first solution of the previous section so it 
describes exactly the same physics. 

4 Conclusion 

We now summarize our two main results, with their implications. Firstly, for each of 
the models, the free energy (up to overall scalings) was the same as other models previ- 
ously studied and using completely different matrix ensembles, as revealed by the ratio 
^2§2. = —13/9) This means the following: if these models were actually describing gravity 
coupled to some other system, then we expect that different regulators would introduce 
a dependence of the free energy on a coupling parameter. Hence using different matrix 
ensembles would yield different results. But this is not the case so our analysis confirms 
that all of these models describe pure gravity including both oriented and unoriented 
surfaces, as it was first assumed for previous solutions. 

Secondly, although both models yield exactly the same result, there is an additional 
solution for n-odd (ratio 83/9) . This solution differs from the other one only by the 
coefficient of the torus/Klein bottle term. So all ratios that do not involve this term are 
the same in both solutions. A similar result was found in p] for QRSD matrices, where 
one solution did not include odd Euler character surfaces. In the present case, the physical 
interpretation of the extra solution remains unclear. 

So far, the C I generator ensemble remains unstudied (as are most of the circular en- 
sembles, which integrate over the entire symmetric space [Q). An interesting observation 
is that measures for D III - n even - and C I are similar to A II and A I respectively after a 
change of variables y = x"^. The D III and C I generator ensembles are more complicated 
to analyse though, because their range of integration is from to 00 instead of —00 to 
00, which causes the appearance of boundary terms in the recursion relations. Given that 
D III with n even and A II yield the same physics in the double-scaling limit, we might 
expect C I and A I to do as well. 
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Cartan's 






System of 




notation 


X = G/K 


Rank 


restricted 
roots 


Multiplicity 


AI 


SL{n,R)/SO{n) 


11- 1 


At-i 


iria = 1 


All 


SU{2n)/Sp{n) 


n- 1 


^n-l 


77l„ = 4 


A III 


SU{p,q)/S{U{p)xU{q)) 


n=mm{p,q) 


P=l,Cn 

p > q,BCn 


7?7„ = 2, 7772/3 = 1 
lUa = 2, 7772/3 = 1 

777/3 = 2(p - q) 


BDI 


S0ip,q)/S0ip)xS0iq) 


n=mm{p,q) 


P = q,Dn 
P > q,Bn 


rua = 1 

mp = {p- q) 

lUa = 1 


D III 


S0{2n)/U{n) 


n/2 


n = 2k,Ck 
n = 2k + l 
BCk 


777c, = 4, 7772/3 = 1 
777c, = 4, 7772/3 = 1 
777/3 = 4 


CI 


Sp{n,R)/U{n) 


n 


Cn 


777,, = 1 


CII 


Sp{p,q)/Spip)xSpiq) 


n=mm{p,q) 


p = q,Cn 
p>q, BCn 


777q, = 4, 7772/3 = 3 
777q, = 4, 7772/3 = 3 

777/3 = 4(p - g) 



Table 1: Classification of the symmetric spaces. 
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Cart an 's 


System of 


Jacobian for 


notation 


restricted roots 


generator ensemble 


AI 


An-l 


111*^* '^j 


All 


An-l 




A III 


P=l,Cn 


m^'-^'ni^^'^'-'^^' 




p> q , BCn 


m-' - -'ru -^\'^'-''^' 


BDI 


p = q , Dn 






p> q , Bn 


':\A-AVi.\^^'~' 


CI 


^n 


11 a;' -xj l[\xi 


CII 


p = q ,Cn 


fiU-x'.rii^^' 




p> q , BCn 


m4 - x^rii x^\'^'-'^'-' 


Dili 


n = 2k ,Ck 


"'n(^.-^?rn^^ 




n = 2k + l, BCk 





Table 2: Symmetric spaces and their Jacobians. 



order 


eq.(D 


eq.® 


/9° 


4/9 = 4/9 




1=1 




p-f" 


9o = 2f 




9o = 2f 




^-M-^ 


91 = f'/2 




9o = 4^1 




p-2^ = /?-/^-2^ 


g, = 1/16/" - p 


4t = 2/3 p - 


- 1/3 92 + 1/12 g[- 


-1/48(7^' 



Table 3: Lowest order solutions of the two recursion relations for R2N+1 and S'2Ar+;. 
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order 


eq-dH) 


eq.(|D 


r'/' 





hoXc = 


^-2/5 


6/X, = 


-6/ - 5fo^c + hXc - hoko - k'^ = 


^-3/5 


+/i'o' - 6/A;o - 6//ioXe = 


3/' + 6fho + goko - hiko + hok'^ - hoki 
k'^/2 - giX, + g'^X, + h^X, -k[ = Q 



Table 4: Lowest order solutions of the two recursion relations for the Dili - odd n 
matrix model. 



order 


eq.(^) 


^-2/5 






Table 5: Lowest order solutions of the recursion relation for the Dili - even n - matrix 
model. 



